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ABSTRACT 

We report on a water maser survey towards a sample of 27 planetary nebulae 
(PNe) using the Robledo de Chavela and Medicina single-dish antennas, as well 
as the Very Large Array (VLA). Two detections have been obtained: the already 
known water maser emission in K 3-35, and a new cluster of masers in IRAS 
17347-3139. This low rate of detections is compatible with the short life-time of 
water molecules in PNe (~ 100 yr). The water maser cluster at IRAS 17347- 
3139 are distributed on a ellipse of size ~ 0'/2 x O'/l, spatially associated with 
compact 1.3 cm continuum emission (simultaneously observed with the VLA). 
From archive VLA continuum data at 4.9, 8.4, and 14.9 GHz, a spectral index 
a = 0.76 ±0.03 (S u oc v a ) is derived for this radio source, which is consistent with 
either a partially optically thick ionized region or with an ionized wind. However, 
the latter scenario can be ruled out on mass-loss considerations, thus indicating 
that this source is probably a young PN. The spatial distribution and the radial 
velocities of the water masers are suggestive of a rotating and expanding maser 
ring, tracing the innermost regions of a torus formed at the end of the AGB 
phase. Given that the 1.3 cm continuum emission peak is located near one of 
the tips of the major axis of the ellipse of masers, we speculate on a possible 
binary nature of IRAS 17347-3139, where the radio continuum emission could 
belong to one of the components and the water masers would be associated with 
a companion. 

1 Laboratorio de Astrofisica Espacial y Fi'sica Fundamental, INTA, Apartado 50727, E-28080 Madrid, 
Spain; e-mail: itziar@laeff.esa.es; jfg@laeff.esa.es 

2 National Radio Astronomy Observatory, PO Box O, Socorro, NM 87801 

3 Centro de Radioastronomfa y Astrofisica, UNAM, Apdo. Postal 3-72 (Xangari), 58089 Morelia, Mi- 
choacan, Mexico 

4 Instituto de Astrofisica de Andaluci'a (CSIC), Apartado 3004, E-18080 Granada, Spain 
5 INAF, Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy 
6 Instituto de Ciencias del Espacio (CSIC)-IEEC, Gran Capita 2, E-08034 Barcelona, Spain 



- 2 - 



Subject headings: planetary nebulae: general — planetary nebulae: individual 
(K 3-35, IRAS 17347-3139) — masers — radio lines: ISM 

1. Introduction 

Nowadays we know that planetary nebulae (PNe) are the result of the evolution of 
carbon- and oxygen-rich red giants, after going through the AGB and a brief (~ 1000 yr) 
proto-PN phase (Kwok 1993). In a simplified picture, a PN forms when the remnant 
(central) star reaches an effective temperature (T e g) of ~ 30 000 K and the envelope ejected 
during the AGB phase is ionized. However, the transformation of an AGB star into a PN 
is much more complex. The spherical symmetry typical of AGB envelopes is lost in the 
transition, and most PNe exhibit elliptical and bipolar shells (Balick 1987). In addition, 
highly collimated outflows develop in the proto-PN and/or PN phase, which may play a 
crucial role in the dynamical and shaping processes of PNe (Sahai & Trauger 1998). 

Envelopes of evolved stars with oxygen-rich chemistry are characterized by strong OH, 
H2O and SiO maser emission (Reid & Moran 1981; Elitzur 1992; Habing 1996). Because of 
its high intensity and the possibility to be observed at high spatial and spectral resolution, 
maser emission is ideal to study dynamical and physical conditions of circumstellar envelopes 
(Habing 1996; Engels 2002). Masers are stratified in these envelopes with SiO masers located 
close to the stellar surface, water masers at about 10-100 AU, and OH masers further away, 
up to ~ 10 4 AU (Lane et al. 1987). In the transformation of an AGB star into a PN, 
masers disappear sequentially (Lewis 1989). When pulsation and mass loss cease at the tip 
of the AGB, SiO masers disappear (Gomez et al. 1990; Nyman, Hall, & Olofsson 1998; 
Engels 2002), while water and OH masers can still be found in the proto-PN phase (Likkel 
& Morris 1988; Marvel & Boboltz 1999; Sahai et al. 1999a,b). As the star enters its PN 
phase, water molecules are rapidly destroyed by the ionizing radiation (Lewis 1989; Gomez 
et al. 1990). Thus, water masers are not expected in PNe, and only OH masers seem to 
persist at the very early stages of their evolution. 

Surprisingly, water maser emission has recently been found in the PN K 3-35 (Miranda 
et al. 2001, hereafter MGAT01). K 3-35 is a young bipolar PN that contains a bright 
compact core, a dense equatorial torus and a bipolar precessing jet, with all these structures 
surrounded by a faint elliptical envelope (Aaquist & Kwok 1989; Aaquist 1993; Miranda 
et al. 1998, 2000, 2001). Water maser emission was detected in two different regions of K 
3-35: in the torus, at a distance of about ~ 85 AU from the center (assuming a distance 
of 5 kpc to K 3-35, Zhang 1995), and at the tips of the bipolar radio jet, at ~ 5,000 AU 
from the central star. The time scale over which water molecules are expected to survive 
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destruction by the ionization front is very short, ~ 100 years (Gomez et al. 1990). Thus, 
it was inferred that K 3-35 had been caught at the very moment of its transformation into 
a PN (MGAT01). In addition, water maser emission was not expected so far away from the 
central star, because of the lack of physical conditions to provide the energy to pump the 
water maser (Marvel 1997). The bipolar jet in K 3-35 was proposed as the excitation agent 
of the distant water masers, because of the former's spatial and kinematic association with 
the latter. For this source, some unknown and local shielding mechanism against ionizing 
radiation was invoked to explain why water molecules still survive. Some recent models of 
gas-phase chemistry, however, predict the production of water and other molecules in the 
envelopes of PN, long after the original molecules were photodissociated (Ali et al. 2001; 
Hasegawa et al. 2000). 

To better understand the phenomenon of water emission in PNe, and to determine 
whether K 3-35 is representative of this kind of sources, or an extremely rare case, we have 
undertaken a survey of water maser emission towards PNe. In this paper we present the 
results of our survey. A technical description of the observations is presented in §2. Our 
results are presented in §3 and their discussion in §4. Finally, we summarize our conclusions 
in §5. 



2. Observations 

We observed the 6\q ^§23 transition of the water molecule (rest frequency 22235.080 
MHz) towards a sample of planetary nebulae (see Table 1), using the NASA 70-m antenna 
(DSS-63) at Robledo de Chavela (Spain), the 32-m antenna 7 at Medicina (Italy), and the 
Very Large Array (VLA) of the National Radio Astronomy Observatory 8 (USA). Several 
sources were observed with more than one telescope. Further simultaneous water maser and 
1.3 cm continuum observations were made with the VLA towards IRAS 17347-3139, the 
source where maser emission was first detected at Robledo de Chavela. We also analyzed 
VLA archive continuum data of this source. Technical details of these observations are 
described below. 

The observed sources were chosen with the criterion of showing some sign of being 
associated with molecular gas and/or dust, e.g., CO or OH lines, or (far or near) infrared 



7 The Medicina VLBI radio telescope is operated by the Radioastronomy Institute of the Italian National 
Research Council (CNR) 

8 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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emission. 



2.1. Robledo de Chavela 70-m 



The 1.3 cm receiver of this antenna was a cooled HEMT. Water maser observations were 
carried out using a 256-channel autocorrelation spectrometer, covering a bandwidth of 10 
MHz, which provides a velocity resolution of 0.5 km s -1 . At this frequency, the half-power 
beamwidth of the telescope is ~ 41". Spectra were taken in frequency switching mode, with 
a switch of 5 MHz, thus providing an effective velocity coverage of 202.4 km s -1 (15 MHz) 
centered at the VLsr of each source. Only left circular polarization was processed. Calibration 
was performed using a noise diode. These observations were carried out between 2002 April 
and August, with system temperatures between 45 and 135 K, and a total integration time 
of 20 min per source. The rms pointing accuracy was better than 15". A total of 15 sources 
were observed with this antenna (see Table 1). The data reduction was performed using 
the CLASS package, which is part of the GAG software package developed at IRAM and 
Observatoire de Grenoble. 



2.2. Medicina 32-m 



Observations with the Medicina 32-m antenna were carried out in several periods from 
March to October 2002, making use of a digital autocorrelator spectrometer with a band- 
width of 8 MHz (108 km s _1 at the frequency of the H 2 6i 6 -5 2 3 line) and a spectral 
resolution of 9.77 kHz. A GaAs FET amplifier was used as a front-end, which allows for 
zenith system temperatures of ~120 K under good weather conditions. During the obser- 
vations the system temperature varied from 130 to 290 K depending on weather conditions 
and elevation. Observations were made in total power mode, with integrations of 5 min on- 
and 5 min off-source. Then the observation was repeated as many times as needed to reach 
typical la rms noise levels of ~0.5 Jy. For each source, all spectra taken on the same day 
were averaged and smoothed to a spectral resolution of 0.263 km s -1 . The antenna gain 
as a function of elevation was determined by total power integrations on DR21 (adopted 
flux density 18.8 Jy) repeated every 1-2 hours. Pointing was checked on strong water maser 
sources and the resulting accuracy is 20". Nine sources were observed with this antenna. 
The data reduction was performed with CLASS. 
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2.3. Very Large Array 

The VLA observations were made in two different periods. As part of our water maser 
survey, on 2002 May 5 we observed 14 planetary nebulae in the A configuration. The in- 
tegration time for each source was ~ 15 minutes. We sampled 63 channels with 97.7 kHz 
resolution (~1.3 km s _1 at the observed frequency), in both right and left circular polariza- 
tions. The source 3C 48 was used as the primary calibrator for flux density, for which we 
adopted a flux density of 1.13 Jy using the latest VLA values (1999.2). Amplitude calibra- 
tion was performed using the model of 3C 48 downloaded from the NRAO web site. Phase 
calibrators used for each target source are listed in Table 2. The calibration and data reduc- 
tion were carried out using the Astronomical Image Processing System (AIPS) of NRAO. 
Maps were produced setting the "robust" weight parameter to 0, to optimize the tradeoff 
between angular resolution and sensitivity, and deconvolved using the CLEAN algorithm. 
The resulting synthesized beams are also shown in Table 2. In K 3-35, self-calibration was 
performed, given its strong water maser emission. To mitigate the Gibbs ringing, spectral 
Hanning smoothing was applied to all sources, giving a final velocity resolution of 2.6 km 
s- 1 . 

On 2002 July 9 we made simultaneous observations of 1.3 cm continuum and water 
maser towards IRAS 17347-3139 with the B configuration, using the 4IF spectral line mode 
of the VLA. For the continuum measurements, we used two of the four IFs, with a bandwidth 
of 25 MHz, and seven channels 3.125 MHz wide each, centered at 22285.080 MHz. For the 
maser emission observations (made with the other two IFs) we used a bandwidth of 3.125 
MHz, sampled with 63 channels of 48.8 kHz each (0.66 km s" 1 ). Spectral Hanning smoothing 
was applied, with a final velocity resolution of 1.3 km s _1 . The source 3C 286 was used as 
the primary calibrator for flux density (adopted flux density = 2.53 Jy), while the source 
J1745-290 was used as phase calibrator (bootstrapped flux density = 0.986±0.008 Jy). The 
on-source integration time was ~ 1 h. Initial calibration was done separately for each pair of 
IFs to remove electronic phase differences between IFs. Self-calibration was performed using 
the strongest water maser component identified in the narrow bandwidth. The phase and 
amplitude corrections were then applied to both the narrow and broad bandwidth, removing 
both atmospheric and instrumental errors. Cleaned maps for both line and continuum were 
produced by setting the "robust" weight parameter to 0, resulting in a synthesized beam of 
0"81 x 0"26 (PA = 26°). 

Finally, we have analyzed VLA archive continuum data at 4.9, 8.4, and 14.9 GHz in IRAS 
17347-3139 (from VLA project AP192). These data were observed in DnC configuration on 
1991 February 24. Two independent IFs of 50 MHz bandwidth were used, providing an 
effective bandwidth of 100 MHz. Both right and left circular polarizations were observed. 
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The on-source integration time at each of the observed frequecies was ~ 10 minutes. The 
source 3C 286 was used as flux calibrator, with adopted flux densities of 3.45 Jy (14.9 
GHz), 5.21 Jy (8.4 GHz), and 7.48 Jy (4.9 GHz). The phase calibrator was B1730-130, 
with bootstrapped flux densities of 7.02±0.20 Jy (14.9 GHz), 6.66±0.06 Jy (8.4 GHz) and 
6.63±0.02 Jy (4.9 GHz). Data were self calibrated, and cleaned images were obtained setting 
the "robust" weight parameter to 0, resulting in synthesized beams of 4'.'24 x 3713 (P.A. = 
44°), 7"90 x 5"56 (P.A.= 40°) and 12"62 x 10755 (P.A.= 61°) at 14.9, 8.4, and 4.9 GHz 
respectively. 

3. Results 

3.1. Survey of water masers 

Table 1 shows the result of our survey of water maser emission towards a sample of 
27 planetary nebulae. Emission was detected in two sources of the sample. One of them, 
K 3-35, was previously known to harbor a water maser, being the first confirmed case of a 
bona-fide planetary nebula showing this kind of emission (MGAT01). We have also obtained 
a new detection towards IRAS 17347-3139 with the Robledo antenna, on day 131 of 2002 
(and confirmed with the VLA, on day 189 of 2002, see Fig. 1), which could represent the 
second case of a planetary nebula associated with water maser emission. 

3.2. IRAS 17347-3139 

3.2.1. Radio continuum emission 

Radio continuum emission associated with IRAS 17347-3139 is detected at all four ob- 
served frequencies, unresolved at 4.9, 8.4, and 14.9 GHz, and marginally resolved at 22.3 GHz 
(deconvolved size ~ 073). This continuum emission is located at a(J2000) = 17 fe 38 m 0O!586, 
5(J2000) = -31°40'55'. / 67 (estimated from the 22.3 GHz image, absolute position error 
~ 0725). From the measured flux densities (Table 3), we derive a spectral index (S u cx u a ) of 
a = 0.76 ± 0.03 between 4.9 and 14.9 GHz (see Fig. 4), indicating a thermal nature for this 
radio source. The measured 22.3 GHz flux density is a factor of two lower than its extrapo- 
lation using this spectral index (97 ± 14 versus ~ 200 mJy). Part of this loss of flux density 
could be due to source variability, since the 22.3 GHz data were observed at a different 
epoch. Another possibility is that the 22.3 GHz observations, being obtained with a more 
extended VLA configuration than at lower frequencies, resolve out weak, extended emission 
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due to the lack of short spacings in the interferometer, thus resulting in a net decrease of 
flux density. We also note that these data were observed at low elevation, which makes 
the absolute flux calibration at this wavelength more uncertain. However, it is also likely 
that a significant fraction of the difference can be explained by the spectral index becoming 
flat between 14.9 and 22.3 GHz. This could indicate that the continuum emission becomes 
optically thin between these frequencies. Quasi-simultaneous, matching-beam observations 
around this frequency range would be needed to test this possibility, ruling out variability, 
resolution, and calibration effects. 

Extended and weaker radio continuum emission, located ~ 45" southeast from the 
unresolved source, is also detected in the 4.9 and 8.4 GHz maps. With such a large separation, 
this extended emission is not likely to be related to IRAS 17347-3139. 

3.2.2. Water maser emission 

Thirteen water maser spots were detected with the VLA on day 189 of 2002 (see Table 
4 and Fig. 2). The spectrum of the whole maser emission, integrated over an area of 
l'/3 x 2'.'3, can be seen in Fig. lb. The maser spots trace an elliptical structure, with the 
radio continuum peak near one of the tips of its major axis. A least-square fit of the maser 
positions to an ellipse gives axes of ~ 0'/25 x 0'/12, with its major axis oriented at PA ~ 62°. 
From the axis ratio of the ellipse, if we assume that this ellipse is tracing a circular ring, the 
axis of the ring would be oriented at ~ 60° with respect to the line of sight. 

In addition, an analysis of the variation of the radial velocities in the ellipse reveals some 
trends. If we assume that the mean value of the observed velocities (—65.1 km s -1 ) represents 
the "systemic velocity" of the maser spots, the maximum (redshifted) radial velocities of ~ 
3-4 km s _1 tend to cluster towards the west of the maser distribution, while the maximum 
(blueshifted) velocities of ~ 3 km s™ 1 are at the east. This would suggest the presence 
of rotational motions, although the kinematic trend is not compatible with pure rotation, 
since the maximum velocities are not exactly at the tips of the major axis of the ellipse. A 
combination of both rotation and expansion in a toroidal structure is necessary to explain 
this kinematic trend (see §4.2.4). 

3.3. K 3-35 

With the VLA, we have detected four maser spots at the central region of this source 
(see Table 4 and Fig 3). In order to compare the positions of the two sets of water maser 
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observations (MGAT01 and this paper), we have obtained a continuum map at 1.3 cm by 
averaging 51 channels free of water maser emission. The peak of the continuum emission 
is located at a(J2000) = 19 /l 27 m 44!025, 5(J2000) = 21°30'03"43 (absolute position error 
~ (X'05), which is consistent, within the errors, with the position reported by MGAT01. The 
two sets of water masers where then aligned with respect to the 1.3 cm continuum peak 
position, assuming it coincident for the two epochs, since the relative positional accuracy 
between continuum and masers within each map (< / /004) is better than their absolute 
position accuracy (~ 0'/05). However, we cannot discard the presence of small variations in 
the ionized structure of the source, that can produce slight changes in the position of the 
continuum peak. The overlay (Fig. 3) shows that the position and velocity distributions 
of the maser spots in 2002 are different from those observed in 1999 by MGAT01, which is 
not surprising giving the typical variability of this kind of emission. In particular, no maser 
emission was detected at the tips of the bipolar lobes of the PN at a 3<r level of 7 mJy, which 
was instead detected in MGAT01. The maser spots observed in 2002 seem to correspond to 
the northwestern masers of the central region C detected in 1999, but shifted by ~ 0'/015 
toward the west. However, the time difference between the two epochs (~ 2.5 yr) is too long 
to make an unambiguous identification of the same water maser spots in the two epochs. 
Therefore a derivation of proper motions is probably meaningless. The observed positional 
variations would also be compatible with destruction of the water maser shell zone observed 
in 1999 by an ionization/ shock front and the creation of a new maser zone ahead of the shock. 
New, more closely spaced multi-epoch observations, and with high S/N for the continuum 
emission (to use it as an accurate position reference for the masers), are crucial to interpret 
the observed positional variations of the water maser spots in K 3-35. 



4. Discussion 

4.1. The lifetime of water molecules in PN 

Two out of 27 PNe observed present water maser emission. This low rate of detection 
seems to indicate that water molecules are destroyed in a very short time after the star 
enters the PN phase, or that the excitation conditions for maser emission cease rapidly. The 
destruction time of water molecules in PNe has been estimated to be ~ 100 yr (Gomez et 
al. 1990). It is difficult to make a quantitative statistical study since we do not know 
the exact age range of our sample (the selection criteria were of qualitative nature and the 
sample is probably biased towards young PNe), and the sensitivities of the telescopes used are 
significantly different. Nevertheless, if we assume that our sample of PNe is homogeneously 
distributed in age, along the ~ 2 x 10 4 yr that the PN phase lasts, our observed detection 
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rate of 1 out of 26 would be compatible with the estimated ~ 100 yr survival time of water 
molecules. 

In the following we will concentrate the discussion on IRAS 17347-3139, whose maser 
emission is reported for the first time in this paper. 

4.2. IRAS 17347-3139 

4.2.1. The distance to IRAS 17347-3139 

As far as we know, there is no published information about the distance to this source. 
A crude estimate for this distance can be obtained taking into account that water maser 
shells in red giants are located at a radius of ~ 10-100 AU (e.g., Spencer et al. 1979; Bowers 
et al. 1993). In addition, there is evidence that the water maser shell moves outwards as 
the star evolves and the mass loss rate increases (Yates & Cohen 1994). Therefore, if the 
angular radius of / /12 of the water maser ring in IRAS 17347-3139 corresponds to a linear 
radius of > 100 AU, a distance of > 0.8 kpc is obtained for the nebula. 

As an independent estimate, from the observed flux density at 4.8 GHz and assuming 
a size of 0'.'3 (the deconvolved size of the 22.3 GHz image) for the ionized region, we can use 
the statistical distance scale for galactic PN proposed by Zhang (1995). We obtain with 
this method a value of ~ 0.8 kpc, which is consistent with the lower limit mentioned above. 

4.2.2. The nature of IRAS 17347-3139 as a planetary nebula 

IRAS 17347-3139 is a source with a large infrared excess (Zijlstra et al. 1989). It 
shows OH maser emission (Zijlstra et al. 1989) but no SiO maser emission (Nyman, Hall, & 
Olofsson 1998). Its IRAS two-color diagram indicates that it is an evolved object (Garcia- 
Lario et al. 1997). However, this diagram cannot discriminate between proto-PNe and PNe. 
The classification of IRAS 17347-3139 as a PN was based on the presence of radio continuum 
emission attributed to an ionized nebula (Zijlstra et al. 1989). Before this classification is 
conclusively established, other possibilities to explain the radio continuum emission should 
be ruled out. 

Strong radio continuum emission at centimeter wavelengths is observed in the Egg Neb- 
ula, a well known proto-PN, although its spectral index (2.6), too large to be due to free-free 
emission, suggests dust as the origin of the radio continuum emission (Jura et al. 2000). 
However, in IRAS 17347-3139, the spectral index of 0.8, being lower than 2, rules out dust 
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as the origin of the radio continuum emission. 

A second possibility to explain the radio continuum emission would be an ionized wind 
from the central star. In order to test this possibility, we estimated the mass-loss rate 
required to explain the observed radio continuum emission as arising from an ionized wind, 
which can be compared with typical values of central stars of proto-PNe and PNe. Assuming 
a spherical wind with a constant velocity of 1000 km s _1 and an electronic temperature of 
10000 K, from the observed flux density at 3.6 cm (Table 4) we obtain a value for the ionized 
mass-loss rate of Mj ~ 1 x 10~ 4 (_D/kpc) 3//2 M yr _1 (Wright & Barlow 1975, Panagia & 
Felli 1975). For the same set of parameters, collimated winds would require smaller mass-loss 
rates. For example, a biconical wind (corresponding to a spectral index of 0.6) would require 
a mass- loss rate ten times smaller, for an angle of aperture of ~ 30°, and a collimated wind 
with a variable angle of aperture (corresponding to the observed value of the spectral index 
of ~ 0.8 between 4.9 and 14.9 GHz) a mass-loss rate three times smaller than a spherical 
wind (see Reynolds 1986). Therefore, in all cases, the derived mass-loss rates are much larger 
than those observed in central stars of PNe and proto-PN (< 10~ 7 M Q yr _1 , Patriarchi & 
Perinotto 1993; Vassiliadis & Wood 1994), which seems to exclude that the radio continuum 
emission arises in an ionized stellar wind. 

All this suggests that the radio continuum emission in IRAS 17347-3139 arises in a 
partially optically thick ionized nebula and that it is probably a very young PN. If we 
assume that the radio continuum emission is optically thin for frequencies above ~ 20 GHz, 
we derive an electronic density of ~ 1 x 10 6 cm -3 for an homogeneous HII region at a distance 
of 0.8 kpc. The ionizing photon rate required to mantain such an HII region is A ~ 10 46 
s _1 , that would correspond to a star with T cff > 26 000 K, according to Panagia (1973). 
This temperature is higher than the T cff ~ 15 000 K deduced by Zijlstra et al. (1989), but 
it is more typical of young PNe. 

Finally, this source presents also infrared lines of [Ne II] and Bra (P. Garcia- Lario 2003, 
private communication), which confirms its identification as a PN. 

4-2.3. Relationship between masers and infrared emission 

It is interesting to compare the water maser data with the K-band image of the object, 
taken with the NICMOS instrument of the Hubble Space Telescope. We have retrieved this 
image from the HST archive and it is shown in Fig. 5. A preliminary description can be 
found in Bobrowsky, Greeley, & Meixner (1999). In addition to a faint envelope of ~ 3'/5 
in size, the brightest regions of this source trace a clear bipolar structure with a size of 
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~ 2'.'0 x 0"7 and the main nebular axis oriented at PA ~ —40°. The NW lobe is brighter 
than the SE one and they are separated by a dark lane that suggests the presence of an 
equatorial torus-like structure. 

We also note that there is a positional shift between the nominal infrared position of 
the nebula given by the HST and the peak position of the radiocontinuum emission, the 
latter being displaced ~ 1" south from the center of the nebula. This angular separation 
might cast some doubts on the relationship of the radiocontinuum and the water masers 
with the infrared nebula. However, the 2MASS catalog lists an infrared point source at 
a(J2000) = 17 /l 38 m 00!610, 5(J2000) = -31°40'55'/23 (absolute position error ~ 0'/06), which 
is ~ 0'.'5 northeast of the radio continuum position, i.e., closer to the latter than the HST 
image center. An offset of ~ 0'.'5 is reasonable even if it is larger than the formal radio 
continuum position error (0'/25), since we do not expect the maximum of the unresolved 
infrared 2MASS image to exactly coincide with the radio emission, taking into account that 
the resolved HST image has its maximum in the northern part of the nebula. Therefore, 
the discrepancy with the HST position is probably due to errors in the HST Guide Star 
Catalog, which can be of the order of 1". Moreover, the association between the radio and 
infrared emissions is reinforced by two further arguments: (1) The thermal nature of the 
radio source (see §4), which seems to exclude the presence of a background extragalactic 
source, and (2) the distribution of the water maser spots, with the major axis of the ellipse 
oriented almost perpendicular to the main nebular axis. This indicates a close relationship 
and that the water maser emission could arise in the dense equatorial ring related to the 
dark lane observed in the K-band infrared image. 



4-2.4- The ring of water masers 

As we mentioned above, to explain the kinematics of the water masers, we need to 
assume the presence of both rotation and expansion in a ring. The uncertainty in the 
distance to the source makes it difficult to estimate the velocity component due to Keplerian 
rotation. If we assume that the expansion and rotation velocity are of the same order, we 
obtain a reasonable fit for a torus radius of 0'/12, with a velocity of 2.4 km s -1 . For such 
a rotational velocity, the central mass should be M* = 0.78D M Q , where D is the distance 
from the Sun to the source, in kpc. This implies a central mass of ~ 0.6 M Q , assuming a 
distance of 0.8 kpc. The water maser ring could trace the innermost regions of the torus-like 
density enhancement formed at the very end of the AGB evolution of the precursor star of 
IRAS 17347-3139. 

As already mentioned, the radio continuum emission peak is located out of the ellipse 
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traced by the water maser spots. An explanation for this could be the presence of inho- 
mogeneities in the ionized region. Nevertheless, we also speculate on the possible nature of 
IRAS 17347-3139 as a binary system, so that the position of the continuum peak represents 
the true central star, while the water maser ellipse is associated with a companion. Our 
observations have not enough resolution to confirm the binary nature of this source, since 
the center of the water maser ellipse and the continuum peak are separated a distance of 
~ C/15, smaller than the beam size (see Fig. 2). This interpretation might be confirmed 
by resolving the ionized region, which would allow us to determine its small-scale geometry, 
and whether a double source is indeed present. 

The water maser emission in both K 3-35 (see MGAT01) and IRAS 17347-3139 (this 
paper) are suggestive of equatorial ring-like structures which have been formed at the end 
of the AGB phase. In the case of IRAS 17347-3139, the kinematics of the ring suggests 
the presence of both rotating and expanding motions, while the observations of K 3-35 are 
compatible with only expansion in the equatorial plane. In addition, the effective tempera- 
ture for the central star of K 3-35 should be > 60 000 K, as indicated by the He II emission 
(Miranda et al. 2000), larger than that obtained for the central star of IRAS 17347-3139 (> 
26 000 K). Therefore, it is possible that IRAS 17347-3139 represents an earlier stage in PN 
formation than K 3-35. The fast wind in IRAS 17347-3139 has not yet been able to sweep 
up enough material in the (dense) equatorial plane and the water maser zone still preserves 
in part its original kinematics. 

5. Conclusions 

In this paper, we present a survey for water maser emission towards a sample of 27 
planetary nebulae, using the antennas at Robledo de Chavela, and Medicina, as well as the 
VLA. Our main conclusions are as follow: 

1. Besides the known water maser in K 3-35 (MGAT01), we have detected a new source 
of maser emission towards IRAS 17347-3139 using the Robledo antenna (and later 
confirmed it with the VLA). 

2. The low rate of new detections (1 out of 26) is compatible with the assumption that 
water molecules are rapidly destroyed (~ 100 yr) after the star enters the PN phase. 

3. Water masers in IRAS 17347-3139, observed with the VLA, show an elliptical distribu- 
tion. This spatial distribution, and their radial velocities are suggestive of a rotating 
and expanding maser ring, tracing the innermost regions of a torus formed at the end 
of the AGB phase. 
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4. The continuum emission in IRAS 17347-3139 seems to arise from a partially optically 
thick ionized region, since the mass loss rate obtained under the assumption of a wind 
is much larger than the values observed in PNe and proto-PNe. 

5. We speculate on the possible nature of IRAS 17347-3139 as a binary star, where the 
peak of the radio continuum emission represents one of the components, and the maser 
ellipse is associated with a companion. 
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Fig. 1. — Water maser spectra of IRAS 17347-3139 (a) obtained with the Robledo antenna 
on day 131 of 2002, and (b) with the VLA (integrated over a region of l'/3 x 2 / /3, which 
includes all the detected maser components) on day 189 of 2002. 
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Fig. 2. — IRAS 17347-3139 water maser positions. Crosses mark the positions of the VLA 
maser components. The filled circle marks the position of the 1.3-cm continuum emission 
peak [a(J2000) = 17 fc 38 m 00!586, 5(J2000) = -31°40'55"67], which is the reference position 
in this map (relative position error between masers and continuum ~ 0'/03). The beam size 
is 0"81 x C/26 (P.A = 26°) 
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Fig. 3. — K 3-35 water maser distribution. Triangles mark the positions of the maser 
components obtained in our new VLA observations, crosses mark the positions obtained 
in MGAT01 (see Gomez et al 2003). The LSR velocity (in km s _1 ) of each component is 
shown. The filled circle marks the position of the 1.3-cm continuum emission peak, assumed 
equal for both sets of observations, and it is the reference position for this map. Relative 
position errors between maser components and continuum emission are ~ 0'.'004 for our new 
observations, and ~ O'.'OOl for those of MGAT01. The beam size is 0716 x 0'.'09 (P.A = 52°). 
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Fig. 4. — Log-log plot of flux density of IRAS 17347-3139 as a function of frequency, including 
2(7 error bars. The solid line is a linear fit between 4.86 and 14.94 GHz, which gives the 
spectral index. 
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Fig. 5. — Infrared image (K band) of IRAS 17347-3139 taken with the Hubble Space Tele- 
scope. Grey levels are arbitrary, and have been chosen to show the main structure of the 
object. Absolute coordinates of the HST image are probably in error by ~ 1". The filled cir- 
cle shows the nominal position of the 22.3 GHz radio continuum (positional uncertainty 0V25. 
The triangle marks the position of the 2MASS infrared point source (positional uncertainty 
0"06) 
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Table 1. Water maser survey. 



Source 


Right Ascension 3 - 


Dcclination a 


V LSR b 


V ■ c 

v mm 


V c 

v max 


Flux dcnsity d 


DoY c 


Telescope 




(J2000) 


(J2000) 


(km s -1 ) 


(km s- 1 ) 


(km s- 1 ) 


(mjy) 


(of 2002) 




PK 130-10.1 


01 42 19.90 


51 34 37.0 


(-21.5) 


-122.7 


+79.7 


<400 


225 


Robledo 


K 3-94 


03 36 08.10 


60 03 46.6 


(-71.0) 


-111.5 


-30.5 


<11 


125 


VLA 


IC 2149 


05 56 23.90 


46 06 17.0 


(-23.6) 


-76.2 


+28.9 


<900 


114 


Medicina 








(-38.0) 


-78.5 


+2.5 


<11 


125 


VLA 


NGC 2440 


07 41 55.40 


-18 12 31.0 


(+52.2) 


0.2 


+104.2 


<2200 


275 


Medicina 








(+39.6) 


-13.3 


+92.6 


<1400 


300 


Medicina 


NGC 4361 


12 24 30.90 


-18 47 05.0 


(+20.0) 


-81.2 


+121.2 


<300 


107 


Robledo 








(+20.0) 


-81.2 


+121.2 


<230 


108 


Robledo 








(+33.6) 


-19.8 


+86.9 


<5000 


178 


Medicina 








(+27.0) 


-20.1 


+74.2 


<1200 


276 


Medicina 








(+19.6) 


-34.0 


+73.1 


<1400 


300 


Medicina 


M 2-9 


17 05 37.80 


-10 08 32.0 


(+101.4) 


+47.8 


+154.9 


<1400 


81 


Medicina 








(+87.9) 


+34.9 


+140.9 


<1300 


301 


Medicina 


IRAS 17106-3046 


17 13 51.70 


-30 49 40.0 


(0.0) 


-101.2 


+101.2 


<1000 


218 


Robledo 


NGC 6309 


17 14 03.60 


-12 54 37.0 


(-33.8) 


-87.6 


+20.0 


<1700 


109 


Medicina 








(-34.0) 


-87.4 


+19.4 


<11000 


179 


Medicina 








(-40.0) 


-86.9 


+7.0 


<1100 


276 


Medicina 








(-47.1) 


-100.9 


+ 6.7 


<1400 


301 


Medicina 


IRAS 17207-2856 


17 23 55.80 


-28 59 32.0 


(-18.0) 


-119.2 


+83.2 


<1000 


218 


Robledo 


IRAS 17347-3139 


17 38 00.70 


-31 40 54.0 


-62.1+0.5 


-226.2 


+81.2 


2100+300 


131 


Robledo 








-68.0+0.5 






2100+300 


131 


Robledo 








-68.1+0.5 


-168.2 


+34.2 


700+400 


140 


Robledo 








-61.8+1.3 


-95.7 


-54.3 


1360+60 


189 


VLA 








-68.2+0.5 


-166.2 


+36.2 


1700+600 


218 


Robledo 


IRAS 17375-3000 


17 40 44.80 


-30 02 00.0 


(-25.0) 


-126.2 


+76.2 


<1200 


218 


Robledo 


IRAS 17423-1755 


17 45 14.20 


-17 56 47.0 


(+35.0) 


-66.2 


+136.2 


<190 


108 


Robledo 








(+35.0) 


-66.2 


+136.2 


<800 


218 


Robledo 


NGC 6572 


18 12 06.28 


06 51 12.6 


(-10.0) 


-111.2 


+91.2 


<110 


110 


Robledo 








(-10.0) 


-50.5 


+30.5 


<12 


125 


VLA 


Cn 3-1 


18 17 34.09 


10 09 03.4 


(+38.0) 


-63.2 


139.2 


<140 


110 


Robledo 








(+38.0) 


-2.5 


+78.5 


<13 


125 


VLA 


M 3-27 


18 27 48.27 


14 29 06.2 


(-24.0) 


-64.5 


+ 16.5 


<13 


125 


VLA 


M 2-43 


18 26 40.10 


-02 42 57.0 


(0.0) 


-101.2 


+101.2 


<140 


108 


Robledo 








(0.0) 


-101.2 


+101.2 


<300 


167 


Robledo 








(+11.8) 


-40.9 


+64.4 


<1200 


274 


Medicina 


M 1-57 


18 40 20.30 


-10 39 47.0 


(+92.0) 


-9.2 


193.2 


<110 


108 


Robledo 








(+92.0) 


-9.2 


193.2 


<130 


131 


Robledo 


Hu 2-1 


18 49 47.56 


20 50 39.4 


(+46.8) 


-6.5 


+100.0 


<900 


81 


Medicina 








(+33.0) 


-7.5 


+73.5 


<13 


125 


VLA 


IC 4846 


19 16 28.23 


-09 02 36.6 


(-168.0) 


-269.2 


-66.8 


<100 


108 


Robledo 








(-168.0) 


-269.2 


-66.8 


<130 


131 


Robledo 








(-156.2) 


-208.8 


-103.5 


<1700 


274 


Medicina 


Vy 2-2 


19 24 22.20 


09 53 56.0 


(-64.0) 


-104.5 


-23.5 


<14 


125 


VLA 


K 3-35 


19 27 44.02 


21 30 03.4 


22.4+0.5 


-76.2 


+126.2 


270+50 


108 


Robledo 








22.6+1.3 


-20.5 


+60.5 


933+5 


125 


VLA 








22.4+0.5 


-76.2 


+126.2 


1330+90 


131 


Robledo 








22.4+0.5 


-76.2 


+126.2 


1410+130 


140 


Robledo 








22.4+0.5 


-76.2 


+126.2 


1330+120 


246 


Robledo 
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Table 1 — Continued 



Source 


Right Asccnsion a 
(J2000) 


Dcclination a 
(J2000) 


V LS R b 
(km s- 1 ) 


V ■ c 
(km s- 1 ) 


V c 
y max 

(km s- 1 ) 


Flux density d 
(mjy) 


DoY e 
(of 2002) 


Telescope 








22.4+0.5 


-76.2 


+126.2 


2100+400 


284 


Robledo 


M 2-48 


19 50 28.00 


25 54 28.0 


(o.o) 


-40.5 


+40.5 


<13 


125 


VLA 


NGC 6884 


20 10 23.59 


46 27 39.4 


(-21-0) 


-61.5 


+19.5 


<12 


125 


VLA 


NGC 6881 


20 10 52.40 


37 24 42.0 


(-14.0) 


-115.2 


+87.2 


<150 


110 


Robledo 








(14.0) 


-54.5 


+26.5 


<13 


125 


VLA 


IC 4997 


20 20 08.80 


16 43 53.0 


(-66.0) 


-106.5 


-25.5 


<15 


125 


VLA 


IC 5117 


21 32 30.90 


44 35 47.0 


(0.0) 


-101.2 


+101.2 


<200 


131 


Robledo 








(0.0) 


-40.5 


+40.5 


<12 


125 


VLA 


KjPn 8 


23 24 10.45 


60 57 30.5 


(-21.8) 


-75.4 


+31.8 


<1000 


109 


Medicina 








(-35.0) 


-75.5 


+5.5 


<12 


125 


VLA 



a Coordinatcs of phase center in VLA observations, or pointing position of single-dish data. Units of right ascension are hours, 
minutes, and seconds. Units of declination are degrees, arcminutes, and arcseconds 

b Velocity of the strongest detected maser components. For non-detections, the central velocity of the observed bandwidth is given 
between parentheses 

c Velocity range covered by the observational bandwidth 

d For non-detections, upper limits are 3o\ For detections, uncertainties are 2cr 
e Obscrvation date, day of year in 2002 
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Table 2. Phase calibrators and synthesized beams of VLA water maser observations 3. 



Source 


Phase calibrator 


Flux density b 

(Jy) 


Beam size c 


Beam P.A. C 


K 3-94 


J0244+624 


1.23±0.05 


O'/IO 


X 


0707 


36° 


IC 2149 


J0555+398 


2.62±0.12 


O'/ll 


X 


070 7 


57° 


IRAS 17347-3139 


J1745-290 


0.986±0.008 


0'.'81 


X 


0726 


26° 


NGC 6572 


J1849+005 


0.723±0.003 


0'/16 


X 


0709 


45° 


Cn 3-1 


J1849+005 


0.723±0.003 


Of! 16 


X 


0709 


48° 


M 3-27 


J1849+005 


0.723±0.003 


0'.'16 


X 


0709 


53° 


Hu 2-1 


J1925+211 


2.979±0.130 


0'.'15 


X 


0709 


56° 


Vy 2-2 


J1950+081 


0.334±0.014 


cy/16 


X 


0709 


48° 


K 3-35 


J1925+211 


2.979±0.130 


0716 


X 


0709 


52° 


M 2-48 


J2023+318 


0.634±0.025 


0713 


X 


0709 


59° 


NGC 6884 


J2012+464 


0.485±0.018 


0711 


X 


0709 


-85° 


NGC 6881 


J2015+371 


3.56±0.11 


0713 


X 


0709 


73° 


IC 4997 


J2031+123 


0.529±0.017 


0718 


X 


0710 


56° 


IC 5117 


J2202+422 


2.33±0.08 


0711 


X 


0709 


90° 


KjPn 8 


J2322+509 


0.679±0.020 


0711 


X 


0708 


-45° 



a All sources were observed with A configuration, except IRAS 17347-3139 that 
was observed with B configuration 

b Bootstrapped flux density of the phase calibrator 

c Size and position angle of the synthesized beams of the images 
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Table 3. IRAS 17347-3139 continuum emission 



Frequency 


Flux density a 


(GHz) 


(mJy) 


4.86 


61.67±1.2 


8.44 


93.9±1.8 


14.94 


145±4 


22.29 


97±14 



a Uncertainties are 2a 
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Table 4. Water maser components detected with VLA 



Source 


V LS R a 
(km s- 1 ) 


Flux density 15 
(mjy) 


Right Ascension 
(J2000) 


Declination 
(J2000) 


Position uncertainty* 3 ' c 
(") 


IRAS 17347-3139 


-61.2 


440±50 


17 38 00.589 


-31 40 55.66 


0.04 




-61.8 


1310±60 


17 38 00.5903 


-31 40 55.632 


0.017 




-62.5 


1360±60 


17 38 00.5907 


-31 40 55.628 


d 




-63.1 


720±50 


17 38 00.594 


-31 40 55.62 


0.03 




-63.8 


540±50 


17 38 00.598 


-31 40 55.60 


0.04 




-64.5 


410±50 


17 38 00.601 


-31 40 55.58 


0.05 




-65.1 


260±50 


17 38 00.603 


-31 40 55.57 


0.07 




-65.8 


180±50 


17 38 00.597 


-31 40 55.62 


0.10 




-66.4 


150±50 


17 38 00.589 


-31 40 55.71 


0.13 




-67.1 


330±50 


17 38 00.603 


-31 40 55.64 


0.06 




-67.8 


1020±50 


17 38 00.6052 


-31 40 55.615 


0.021 




-68.4 


940±50 


17 38 00.6041 


-31 40 55.626 


0.022 




-69.1 


220±50 


17 38 00.597 


-31 40 55.73 


0.09 


K 3-35 


23.9 


351±6 


19 27 44.02296 


21 30 03.4404 


0.0012 




22.6 


933±5 


19 27 44.02308 


21 30 03.4412 


e 




21.3 


714±6 


19 27 44.02322 


21 30 03.4422 


0.0006 




20.0 


125±6 


19 27 44.0234 


21 30 03.444 


0.004 



a Velocity of water maser emission 
b Uncertainties are 2a 

c Relative position uncertainties with respect to the reference feature used for self-calibration 
d Reference feature. Absolute position error is ~ 0"25 
e Reference feature. Absolute position error is ~ 0^05 



